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ABSTRACT 

We combine IR, optical and X-ray data from the overlapping, 9.3 square degree NOAO Deep Wide- 
Field Survey (NDWFS), AGN and Galaxy Evolution Survey (AGES), and XBootes Survey to measure 
the X-ray evolution of 6146 normal galaxies as a function of absolute optical luminosity, rcdshift, and 
spectral type over the largely unexplored redshift range 0.1 55 z 0.5. Because only the closest or 
brightest of the galaxies are individually detected in X-rays, we use a stacking analysis to determine 
the mean properties of the sample. Our results suggest that X-ray emission from spectroscopically late- 
type galaxies is dominated by star formation, while that from early-type galaxies is dominated by a 
combination of hot gas and AGN emission. We find that the mean star formation and supermassive 
black hole accretion rate densities evolve like ~ (1 + z) , in agreement with the trends found for 
samples of bright, individually detectable starburst galaxies and AGN. Our work also corroborates the 
results of many previous stacking analyses of faint source populations, with improved statistics. 

Subject headings: cosmology: observations — galaxy evolution 



1. introduction 

There are three primary sources of galactic X-ray emis- 
sion: diffuse, hot gas, accreting stellar remnants, such as 
X-ray binaries, and accreting supermassive black holes, 
i.e., active galactic nuclei (AGN). Apart from nearby el- 
lipticals, which tend to be dominated by hot gas emission 
(e.g., Forman, Jones, & Tucker 1994), local surveys indi- 
cate that X-ray binaries dominate the flux from "normal" 
galaxies with quiescent nuclei (e.g., Fabbiano & White 
2003, Smith & Wilson 2003, Muno et al. 2004). There are 
two distinct populations of X-ray binaries: the short-lived 
(<; 10 6 yrs), high-mass X-ray binaries (HMXBs) and the 
long-lived (i> 10 8 yrs), low-mass X-ray binaries (LMXBs). 
Because of their short lifetimes, one expects the X-ray 
emission from HMXBs to track the current star formation 
rate (SFR). In contrast, the X-ray emission from longer- 
lived LMXBs should track the integrated stellar mass with 
some time lag due to stellar and binary evolution (e.g., 
Ghosh & White 2001, Ptak et al. 2001, Grimm et al. 
2002). 

These expectations are largely borne out by observa- 
tions. There are strong correlations between the hard 
(> 2 keV) X-ray luminosity (presumably from HMXBs) 
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and other star-formation indicators, such as radio, far-IR, 
B-band, and UV flux for both local (e.g., Fabbiano 1989, 
David, Jones, & Forman 1992) and higher redshift (Seib- 
ert et al. 2002, Bauer et al. 2002) galaxy samples. Many 
studies have found a linear correlation between X-ray lu- 
minosity and the SFR (e.g., Ranalli et al. 2002, Grimm 
et al. 2003, Gilfanov et al. 2004, Persic et al. 2004, and 
Colbert et al. 2004). In particular, Grimm et al. (2003) 
found that the hard (2-10 keV) HMXB X-ray luminosity 
scales with the star formation rate according to 

« - 0.67 x 10 40 (x^V) ^ s_1 - « 

In contrast, the (total band: 0.5-8 keV) X-ray emission 
from LMXBs is well-correlated with the if -band flux of 
galaxies, which in turn is a good tracer of stellar mass. 
Kim and Fabbiano (2004) find that 

L LMXB „ 1Q 40 (^£j ergs B -l | (2) 

where L K * = v K L VK ^ = 2.6 x 10 43 ergs s _1 corresponds 
to M K * = -23.4 mag (Kochanek et al. 2001). 

Because the stellar mass density evolves slowly between 
z = and z ~ 0.5 (e.g., Bell 2004 and references therein), 
we expect little change in the number of LMXBs and as- 
sume that Eqn. @ holds out to z ~ 0.5. On the other 
hand, the star formation rate density rises rapidly with 
look back time, (e.g., oc (1 + z ) 2 - 7±0 - 7 Hogg 2001), so we 
should see a correspondingly rapid evolution in the X-ray 
flux of star-forming galaxies due to the increasing number 
of HMXBs. 

The accretion luminosity density of AGN also rises rapidly 
with redshift (e.g., oc (1 + z) 3 - 2±0 ' 8 , Barger et al. 2005) 
with higher densities of more luminous sources at higher 
redshifts (e.g., Barger et al. 2001, Cowie et al. 2003, Ueda 
et al. 2003, Miyaji 2004, Hasinger 2004, Hasinger, Miyaji, 
& Schmidt 2005). The observed trends suggest that faint 
AGN with characteristic X-ray luminosities approaching 
those of bright starbursts, L x < 10 41 ergs s -1 , may be 
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most abundant at z 1 and, if so, they could produce a 
substantial fraction of the "normal" galaxy X-ray flux at 
these redshifts. 

Given their similar evolution, spectral shapes (e.g., Ptak 
et al. 1999), and potentially similar luminosities, it may be 
difficult to disentangle the X-ray emission from AGN and 
star formation at z <, I. Probing the X-ray evolution of 
normal galaxies between z ~ 0.1 and z ~ 1 is also difficult 
because it requires relatively deep observations over rela- 
tively wide areas. In fact, most of our knowledge about the 
X-ray properties of galaxies still comes from either large- 
area, local surveys or high-redshift, small- volume deep fields. 

This is the second in a series of papers in which we 
attempt to bridge this gap by examining the X-ray evolu- 
tion of galaxies within the 9.3 square degree Bootes field 
of the NO AO Deep Wide-Field Survey (NDWFS; Jannuzi 
& Dey 1999, Jannuzi et al., in preparation, Dey et al. in 
preparation). The XBootes survey (Murray et al. 2005) 
obtained a 5 ksec Chandra mosaic of the entire 9.3 square 
degree field based on 126 ACIS-I pointings. These data 
are too shallow, however, to detect typical galaxies even 
at modest redshifts. Fortunately, the wide area of the sur- 
vey allows us to measure the mean properties of a large, 
representative population of galaxies by "stacking" (aver- 
aging) their X-ray emission (Brandt et al. 2001; Nandra 
et al. 2002; Hornschemeier et al. 2002, Georgakakis et al. 
2003, Lchmer et al. 2005, Laird et al. 2005, Laird et al. 
2006, Lehmer et al. 2007). 

In the first paper in the series, Brand et al. (2005; here- 
after B05) employed the stacking technique to examine 
the X-ray luminosity evolution of w 3300 optically lu- 
minous (~ L*), red galaxies with photometric redshifts 
0.3 < z < 0.9. By constraining the sample to have the 
same evolution-corrected, absolute i?-band magnitude dis- 
tribution at all redshifts (M R < -21.3, (Mr) ~ -22.0), 
it was possible to follow the nuclear accretion histories of 
a consistent population of galaxies throughout this epoch. 
Because these massive, early-type galaxies are the typical 
hosts of powerful quasars at higher redshifts (McLeod & 
McLeod 2001; Dunlop et al. 2003), tracking the decay of 
their nuclear activity to lower redshifts is essential to our 
understanding of the decline in supermassive black hole 
(SMBH) accretion from its peak during the quasar phase 
(z ;> 2) to the present. The observed variation in nuclear 
accretion rate of the > L* red galaxies, L x cx (1 + z) 4±2A , 
is broadly consistent with the behavior of bright, individ- 
ually detected AGN Barger et al. (2001; hereafter BaOl) 
and Barger et al. (2005; hereafter Ba05). 

In the present paper we combine the NDWFS and XBootes 
surveys with redshifts from the AGN and Galaxy Evo- 
lution Survey (AGES, Kochanek et al. 2007) to study 
the X-ray properties of a complete, flux-limited sample of 
galaxies as a function of absolute optical luminosity, red- 
shift, stellar mass, and spectral type. While our general 
procedures are similar to those of other recent stacking 
analyses (e.g., Brandt et al. 2001; Nandra et al. 2002; 
Hornschemeier et al. 2002; Georgakakis et al. 2003, Laird 
et al. 2005, Laird et al. 2006, Lehmer et al. 2007), the 
number of objects we consider (« 6500) is 10-100 times 
larger. The size of our sample allows us to convert the 
short mean exposure time of the XBootes survey (5 ksec) 
into stacked, effective exposure times that are large enough 



to detect galaxies at intermediate redshift (Msec). 

In §2 we describe the data and illustrate our ability to 
detect X-ray emission from the AGES galaxies through the 
stacking approach. In §3 we provide further details of our 
stacking analysis and use Monte Carlo simulations to test 
it and to refine our choices of signal and background aper- 
tures. In §4 we present measurements of the mean evolu- 
tion of all the AGES galaxies. We then consider the radial 
emission profiles of a nearby (z < 0.1) subsample in order 
to test for weak/obscured AGN. We also examine the hard- 
ness and X-ray to optical luminosity ratios as a function 
of redshift, stellar mass, and spectroscopic type (late- type, 
early-type, or AGN). Based on these tests, we evaluate the 
relative contributions of hot gas, LMXBs, HMXBs, and 
AGN emission to the observed signal. We then discuss our 
results and compare them to previous studies. Lastly, we 
determine the evolution of the star formation rate per unit 
galactic stellar mass and the star formation rate density 
based on emission from the late-type galaxy sample and 
trace the evolution of the nuclear accretion rate and accre- 
tion rate density of AGN in the early-type galaxy sample. 
We summarize our findings in §5. Throughout the pa- 
per, we assume a flat ACDM cosmology with f2 mj o = 0.3, 
A ,o = 0.7 and H = 70 km s" 1 Mpc" 1 . 

2. THE X-RAY DATA, THE GALAXY SAMPLE AND 
METHODS 

The NDWFS obtained B w , I, R, and (partial) X-band 
observations of two roughly 9.3 square degree regions. The 
Northern Bootes field has also been surveyed at radio (VLA 
FIRST, Becker et al. 1995 and WSRT, dc Vrics et al. 
2002), far-IR (Spitzer/MIPS, Soifer et al., 2004), mid-IR 
(Spitzer/IRAC, The IRAC Shallow Survey, Eisenhardt et 
al., 2004), near-IR (NDWFS and the Flamingos Extra- 
galactic Survey; Elston et al. 2006), UV (GALEX; Hoopes 
2004), and X-ray (XBootes, Murray et al. 2005) wave- 
lengths. In particular, the XBootes survey covered the 
entire Northern Bootes field with 126 17' x 17' ACIS-I 
images of mean exposure time 5 ksec using the Chandra 
X-ray Observatory (Murray et al. 2005). The survey de- 
tected approximately 3300 point sources and 43 extended 
sources in the field with > 4 counts (4767 sources with > 2 
counts, Renter et al. 2005). 

2.1. X-ray Data 

For our current analysis we use the positions of the 
420,000 photons detected in the XBootes observations, 
with weights correcting for off-axis vignetting and vari- 
ations in exposure time, their estimated energies e ± Se, 
and their positions relative to the pointing center of the 
observations. We limited the analysis to photons with 
0.5 < e < 7 keV, dividing them into total (0.5-7 keV), soft 
(0.5 — 2 keV) and hard (2 — 7 keV) energy bands. The mean 
total, hard, and soft backgrounds measured during these 
observations were approximately 2.8, 1.9 and 0.96 x 10~ 3 
counts arcsec -2 , respectively. 

For the normal galaxies, the average ratio between the 
hard and soft X-ray counts is 0.73. To fit this count 
ratio, we assume an absorbed power-law dN(E)/dE oc 
E~ T (counts s _1 keV -1 ) for the intrinsic X-ray spectrum. 
We then use Portable, Interactive Multi-Mission Simu- 
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lator (PIMMS v3.9dfl to find the best-fit photon index, 
r = 1.29, assuming a fairly low absorption column density 
(N(H) = 4 x 10 20 cm~ 2 ), half of which is the Galactic ab- 
sorption for the Bootes field from Stark et al. (1992). Fac- 
tor of two changes in the assumed column densities change 
r by only a few percent. The spectrum of the "normal" 
galaxies is very similar to the that of the Cosmic X-ray 
Background, Tcxb —1.4 (e.g., Tozzi et al. 2001; Nandra 
et al. 2004) and of resolved populations in the Chandra 
Deep Fields (e.g., Civano, Comastri, & Brusa 2005). Pho- 
ton indices in this range are also typical of both X-ray 
binaries and obscured AGN (1 < T <J 2, e.g., Ptak et 
al. 1999, Muno et al. 2004). Based on our model of 
the composite AGES galaxy X-ray spectrum, we estimate 
that one soft, hard or total band count corresponds to 
1.05 x 10 -15 ergs cm~ 2 s~ 4 , 3.65 x 10~ 15 ergs cm -2 s _1 
and 2.14x 10~ 15 ergs cm~ 2 s _1 , respectively, for an on-axis 
source with a 5 ksec exposure time. 

2.2. The Galaxy Sample 

For our stacking analysis we cross-referenced the XBootes 
catalog of X-ray photon positions (i.e., the event list) with 
those of the galaxies observed by the AGN and Galaxy 
Evolution Survey (AGES, Kochanek et al. 2007). We fo- 
cus on the approximately 6500 galaxies in the 2004 AGES 
main galaxy sample. This sample consists of all galax- 
ies with R < 19.2 mag and a randomly selected 20% of 
galaxies with 19.2 < R < 20 mag. We use Vega magni- 
tudes throughout the paper. The spectroscopic redshifts of 
the galaxies were measured using Hectospec, a robotic 300 
fiber spectrograph for the MMT (Fabricant et al. 1998, 
Roll et al. 1998, Fabricant et al. 2005) and are 95% com- 
plete for R < 19.2; our analysis is not significantly affected 
by the modest level of incompleteness. Although we have 
redshifts for roughly half of the remaining 80% of galaxies 
with 19.2 < R < 20 mag, we do not include them in this 
analysis because they were not targeted based on a simple 
R-band flux-limited criterion. 

Redshifts were measured using two independent data 
reduction pipelines (the Hectospec pipeline at the CfA and 
a modified SDSS pipeline at Steward Observatory), and 
the results were verified by visually inspecting the spectra 
for features consistent with the pipeline redshift. Both the 
pipelines and the visual inspection flagged galaxies with 
spectroscopic signatures of AGN, both obvious broad-lined 
systems and those with emission line ratios characteristic 
of Type II optical AGN. We also carried out a principal 
component analysis of the optical spectra (e.g., Formiggini 
& Brosch 2004, Ferreras et al. 2006). In our analysis, the 
level of star formation was largely characterized by the 
ratio of the second component relative to the first, so we 
divided the sample into early and late-type galaxies at the 
median of the distribution of that ratio for the full AGES 
galaxy sample. 

Rest-frame absolute i?-band magnitudes, Mr, were de- 
termined by combining the redshift measurements with 
the standard ACDM cosmology and the kcorrect algo- 
rithm (Blanton et al. 2003) adapted to the AGES data 
and the NDWFS Bw, I, and R photometric filters. Fig- 
ure [1] shows the distribution of the galaxies as a function 
of redshift and absolute magnitude. The median redshift 

9 http:/ /heasarc. nasa.gov/Tools/w3pimms. html 



of the galaxies is (z) ~ 0.25 with an appreciable tail ex- 
tending to z w 0.6. The change in the density of points 
at fainter magnitudes is due to the random 20% sampling 
of the fainter galaxies. To compensate for this effect, we 
weight the emission from these galaxies by a factor of 5. 

We also determined the rest-frame if -band luminosities, 
Lk, with kcorrect. We use the if -band luminosities to 
estimate the contribution of LMXBs to the observed X-ray 
emission (Eqn. \2§ , and to estimate the stellar mass 

M, ~ (8.0 ±2.0) x 10 10 (-p^- ] M Q , (3) 

based on the results of Colbert et al. (2004, also see Gil- 
fanov 2004) . To explore how late- type galaxies of different 
masses contribute to star formation, we then combine our 
estimates of the stellar mass and the star formation rate 
(determined using Eqn. [T]) to compute the Specific Star 
Formation Rate SSFR = SFR/M* (e.g., Cowie et al. 1996, 
Guzman et al. 1997, Brinchmann & Ellis 2000, Juneau et 
al. 2004, Feulner et al. 2004, Perez- Gonzalez et al. 2005). 

If, instead of star formation, the X-ray emission is pri- 
marily due to accretion onto AGN, then we can estimate 
the SMBH growth rates, Mbh, using the scaling of BaOl, 

M BH = 1-76 x 10- 6 (^) ~* L^N M yr-\ (4) 

where is the AGN luminosity in units of 10 40 ergs s _1 

and e is the efficiency of energy conversion in AGN. Cowie 
et al. (2003) and Ba05 showed that the X-ray luminosi- 
ties of both obscured and unobscured AGN (with 10 42 erg 
s _1 < L x < 10 44 erg s _1 ) are consistent with a radiative 
efficiency of e ~ 0.1 all the way to z — 0. Moreover, the 
results of AGN models do not agree with observations if 
one assumes that lower radiative efficiency AGN dominate 
black hole growth, (e.g., Shankar et al. 2004, Hopkins et 
al. 2007). We will therefore adopt e = 0.1 throughout the 
paper. 

To fit the X-ray luminosities of the AGES galaxies as a 
function of K-band luminosity and redshift, we assume 

L x cxi«(l + z)' 3 . (5) 

In order to estimate the mean star formation density p* 
and nuclear accretion rate density pbh we must extrapo- 
late our results from the magnitude-limited AGES galaxy 
samples to fainter fluxes. By assuming our power-law 
model (Eqn. [5]) holds for fainter sources and combining the 
K-band luminosity functions for late- and early-type galax- 
ies measured by Kochanek et al. (2001), with a change of 
variables from Lk to L x , we compute the X-ray luminos- 
ity density as a function of redshift. The X-ray luminosity 
density can then be converted to star formation or nuclear 
accretion rate densities using Eqns. (1) or (4), respectively. 

We restrict our analysis to a sample of 6146 "normal" 
galaxies with no optical or X-ray evidence of AGN activity. 
When stacked, these galaxies produce a total of 486/669 
hard/soft counts above the background level. In Table 1, 
we compare the emission from the normal galaxies to two 
classes of sources we exclude from the analysis. We elim- 
inated 47 spectroscopically identified AGN (3 broad line 
and 44 narrow line) that were not identified in the XBootes 
Survey (Renter et al. 2005), but these produced a stacked 
flux of only 32/27 hard/soft counts. Much more impor- 
tantly, we excluded 58 galaxies that were identified both by 
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Fig. 1. — The absolute R-band magnitudes of the AGES main 
sample of galaxies as a function of redshift. The transition from 
higher to lower density background points reflects differences in sam- 
pling: all R < 19.2 galaxies are present but only a random 20% of 
galaxies with 19.2 < R < 20 have been sampled. Four of the eight 
absolute magnitude strips of galaxies we analyze are outlined. 



XBootes as X-ray sources and by AGES as spectroscopic 
broad line (24) or narrow line (34) AGN. Although these 
58 AGN represent less than 1% of our sample by number, 
they generate 787/1196 hard/soft counts - nearly twice 
the X-ray flux of the normal galaxies and would therefore 
dominate our results if they were not eliminated from the 
analysis (see Fig. [SJ- Note, however, how similar the com- 
posite X-ray spectrum of the remaining "normal" galaxies 
(rNormai — 1-3) is to that of these AGN (Fagn — 1-3) 
and to that of typical AGN, which account for the bulk of 
Cosmic X-ray Background (CXB; Tcxb — 1-4, e.g., Tozzi 
et al. 2001; Nandra et al. 2004). The remaining "normal" 
galaxies therefore provide a representative sampling of the 
CXB source population, and, in fact, comprise ~ 8.8% of 
the total CXB emission, based on the soft and hard band 
CXB flux estimates of De Luca k Molendi (2004). 

3. METHOD, MONTE CARLO SIMULATIONS AND TESTS 

The basic principle of a stacking analysis is simple. The 
X-ray images of a large number of sources are stacked to 
determine the mean properties of a group of objects that 
are individually undetectable (Brandt et al. 2001; Nandra 
et al. 2002; Hornschemeier et al. 2002, Georgakakis et al. 
2003, Lehmer et al. 2005, Laird et al. 2005, Laird et al. 
2006, Lehme et al. 2007). Quantitatively, if stacking a 
sample of n objects yields N s counts in a signal aperture 
of area A s and Nf, counts in a background annulus of area 
Ah, then the mean number of source counts per object in 
the signal region is 




-40-20 20 40-40-20 20 40-40-20 20 40 
RA (kpc) RA (kpc) RA (kpc) 

Fig. 2. — Total band X-ray images of the stacked early-type, 
late-type, and combined samples of galaxies have been Gaussian- 
smoothed on a scale of 5 kpc and divided into three redshift bins: 
< z < 0.15, 0.15 < z < 0.3, and 0.3 < z < 0.6. We have in- 
cluded only galaxies (extended sources) that had no obvious optical 
signatures of AGN activity. In each panel, the outer circle shows 
the fixed 30 kpc aperture we use to extract the signal and the inner 
circle shows the size of the 50% enclosed energy radius averaged 
over the field and converted to a physical scale based on the mean 
redshift of each bin (fpsF ~ 2, 10, and 15 kpc, respectively). 



The uncertainties in (N) src can be computed either by 
assuming a Poisson distribution of X-ray counts per object, 



1 



1/2 



(7) 



(N) s 



N s - ^N b 



(6) 



or by bootstrap resampling of both the galaxy and X-ray 
photon catalogs (event lists). Adopting the latter method, 
we generate 100 bootstrap samples by randomly drawing- 
new lists of galaxies and X-ray photons (with replacement) 
from the input catalogs until each bootstrap sample con- 
tains the same number of galaxies and X-ray photons as 
the real catalogs. Each mock data set is then subjected to 
the same analysis as the real data. Our error estimates are 
defined by the range encompassing 68% of the bootstrap 
results about the true mean. This approach will produce 
more realistic error bars than the Poisson statistics of the 
stacked image, particularly if the net flux is dominated by 
a small fraction of the objects included in the analysis. 

As a general proof of principle for the subsequent anal- 
yses, Fig. [5] shows images of the stacked early-type, late- 
type and combined samples of galaxies in three redshift 
bins, excluding direct X-ray detections and galaxies with 
optical signatures of AGN activity. Based on the Monte 
Carlo simulations we discuss in the next section, we adopt 
a 30 kpc signal aperture and a 30 to 60 kpc background an- 
nulus about the center of each galaxy. Each panel shows 
this signal aperture and the field-averaged point spread 
function (PSF), i.e., the 50% enclosed energy radius aver- 
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Fig. 3. — A comparison of the mean X-ray luminosities of the 
AGES normal galaxies and of random locations in the XBootes field. 
The mean X-ray luminosities (in units of 10 40 ergs s _1 ) of galaxies 
drawn from bins in the half-integer (integer) absolute magnitude 
strips with Mr < —19 (see Figure [TJ are designated, from faint 
to bright, by open (closed) triangles, squares and pentagons. No 
statistically significant X-ray flux was detected from the intrinsi- 
cally faintest galaxies. Most of the trend here is due to the higher 
redshifts of the more massive (and luminous) galaxies in our flux- 
limited sample. The dashed line shows our expected sensitivity limit 
given the 5 ksec exposure time and the typical number of objects 
in each bin (see §3). Note that through stacking we are able to 
detect average X-ray signals that are much less than the luminosity 
equivalent of one total band count (solid line). 



aged over the field and converted to a physical scale based 
on the mean redshift of each bin. We clearly detect a signal 
from these "normal" galaxies. 

To understand the X-ray evolution of the galaxies quan- 
titatively, we bin them by absolute magnitude, spectro- 
scopic type and redshift. We analyze the galaxies in stag- 
gered, 1 mag wide absolute magnitude strips (see Fig.[T]), 
beginning at integer and half-integer absolute magnitudes 
(i.e., -17.5 > M R > -18.5, -18 > M R > -19, -18.5 > 
M R > -19.5, ... - 20 > Mr > -21). We divide 
each strip into into two redshift bins that contain roughly 
equal numbers of galaxies after correcting for the sparse 
sampling. We also analyze the combined, late-type, and 
early-type galaxy samples separately. For each bin we then 
determine the mean X-ray and rest-frame if-band lumi- 
nosities of the member galaxies. Because the X-ray and 
optical luminosities of galaxies are approximately linearly 
correlated (e.g., David, Jones, & Forman 1992; Shapley et 
al. 2001), the X-ray to optical luminosity ratio should pro- 
vide a measure of the redshift evolution of normal galaxy 
X-ray emission that is essentially independent of galaxy 
(optical) luminosity and stellar mass (Eqn. [3]) . 



In addition to spectroscopically identified AGN, we ex- 
cluded directly detected X-ray sources that did not satisfy 

L x < 6.6 x 10 42 erg s _1 (luminosity limit), (8) 

which corresponds to two total band counts observed at 
the upper redshift limit for the normal galaxies: z = 0.6. 
The most luminous detected sources at lower redshifts, 
with 10 41 erg s" 1 < L x < 10 42 erg s , are within the 
range possible for star forming galaxies, (e.g., Norman et 
al. 2004). Moreover, our bootstrap uncertainty estimates 
are close to the Poisson limits, which means that the X-ray 
flux cannot be dominated by a small number of luminous 
sources. 

As an additional check on our ability to detect galax- 
ies and to diagnose problems with background subtrac- 
tion or contamination, we constructed a null sample to 
compare with the fluxes of the real sources. To create the 
null sample we kept the redshifts and (optical) luminosi- 
ties of the galaxies fixed but assigned them random posi- 
tions in the field, excluding locations closer than 60 kpc 
to an AGES galaxy. We then analyzed the random cat- 
alog in the same manner as the actual data. In Fig. [31 
we compare the luminosities of the real galaxies and the 
"empty" regions. While the real galaxies generate a sig- 
nificant signal in all but the faintest (Mr i> —18) bins, the 
mean luminosities of the random fields are always consis- 
tent with zero. The uncertainties in the random signal 
closely follow the expected trend of <& m i n 47rc2| (z) , where 
®mm ~ 4 x 10~ 17 ergs s" 1 cm" 2 is the sensitivity of the 
observations divided by the square root of the number of 
galaxies in a typical bin. Stacking the galaxies in our lu- 
minosity and redshift bins generally provides <; 25 times 
the sensitivity associated with the 5 ksec exposure time 
for individual sources. Over the range of redshifts we 
probe, this flux limit corresponds to threshold luminosi- 
ties of w 10 38 ergs s _1 to 10 40 ergs s _1 . 

The X-ray properties of our stacked galaxy samples are 
summarized in Tables 2 and 3. Table 2 provides upper 
limits on the hard and soft band X-ray luminosities of 
the optically faintest AGES galaxies that were either not 
detected in X-rays or had net stacked emission with S/N 
< 1 in all X-ray bands. Table 3 provides the data for the 
detected (S/N > 1) galaxies. 

3.1. Monte Carlo Simulations 

We used Monte Carlo simulations to test our analysis 
methods and to refine our choices of parameters. We 
started by generating a mock galaxy catalog using the 
Brown et al. (2001) i?-band luminosity function, a simple 
model for galaxy evolution, and the AGES main galaxy 
sample selection criteria including the sparse sampling. 
The general properties of the resulting synthetic catalog 
are quite similar to those of the real one. The absolute 
magnitude strips are a case in point. From faint to bright 
magnitudes, the real strips with integral, absolute mag- 
nitude boundaries in Fig. [1] contain 471, 1621, 2526 and 
1496 galaxies including the corrections for sparse sampling, 
while the corresponding strips in the synthetic catalog con- 
tain 614, 1630, 2218 and 1527 galaxies. 

To simulate the X-ray data, we began by distributing 
w 350,000 background photons at random positions in the 
field with surface densities of 1.9 x 10 -3 hard (2-7 keV) 
and 0.9 x 10~ 3 soft (0.5-2 keV) photons arcsec -2 , based on 
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the background estimates for the XBootes survey (Kenter 
et al. 2005). We then assigned each galaxy an extended X- 
ray component and a point-like nucleus. For simplicity, the 
extended emission consists solely of soft photons and the 
AGN component consists solely of hard photons. These 
components have mean luminosities of 

i X 40,soft = b{L/L*){l + zf 

£x40,hard = 5(L/L»)(l+^) 4 , (9) 

respectively, in units of 10 40 ergs s _1 , where L/L» is the 
rest- frame (i?-band) optical luminosity in units of the char- 
acteristic luminosity, L„, at the knee of the Schecter lumi- 
nosity function. The redshift scaling and the correlation 
between the X-ray and optical luminosities are motivated 
both by previous studies (e.g., David, Jones, & Forman 
1992; Shapley et al. 2001, Brand et al. 2005) and by our 
own results (see §4). 

Based on the X-ray spectrum of the AGES galaxies, we 
converted the X-ray fluxes to counts using soft and hard 
band flux limits of 1.05 x 10~ 15 ergs cm" 2 s _1 and 3.65 x 
10 -15 ergs cm -2 s _1 , respectively and K-corrections of (1+ 
z)( T ~ 2 \ with T = 1.29 (see §2.1). The hard emission was 
distributed about the galactic centers using a Gaussian 
model of the CXO PSF with a 50% enclosed energy radius 
of 

r PSF = Of.'5 + 6'!o(-^j , (10) 

where D is the randomly assigned distance of the galaxy 
from the CXO pointing center in arcminutes (see the CXO 
Proposers' GuidJ 1 ""!). We modeled the extended emission 
as a Gaussian of dispersion = 10(L/ X*) 1 / 2 kpc con- 
volved with the PSF to get an overall radial profile of 



extended 



(R) oc exp 



-R 2 



(11) 



.2K + ^f). 

In this model the X-ray surface brightness of the galaxies 
is independent of optical luminosity. 

3.2. Radial Emission Profiles 

We first verified that our analysis correctly extracts the 
X-ray emission profiles by analyzing the data for the roughly 
800 galaxies in the mock catalog with z < 0.1 (mean red- 
shift (z) ~ 0.07) and within 10' of the pointing center. We 
consider only these lower redshift sources because, in spite 
of the CXO's high, on-axis resolution, the field-averaged 
PSF still approaches 2 kpc at z ~ 0.07. Fig. 2] shows the 
extracted profile for the total band. Using a 30 kpc signal 
aperture and a 30-60 kpc background subtraction annu- 
lus, we recover the input profile (and background level) 
to within the la bootstrap uncertainties. We carry out a 
similar analysis on the real data in §4.1. 

At higher redshifts, where we cannot spatially resolve 
the emission, we concentrate on extracting the redshift 
evolution of the X-ray luminosity. This involves a trade- 
off between the signal-to-noise ratio (S /N) and the need for 
aperture corrections. Larger apertures capture more sig- 
nal but introduce extra noise from the background, while 
small apertures lose some of the signal and require more 
substantial aperture corrections. Redshift-depcndent bi- 
ases will be minimized by using a fixed angular aperture for 

10 http:/ /cxc. harvard.edu/proposer/POG/html 
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Fig. 4. — Monte Carlo test of the stacking procedure. The solid 
curve shows the input total band emission profile in our Monte Carlo 
simulation, and the open circles show the profile extracted by the 
stacking analysis. The extracted profile agrees with the input profile 
in all radial bins. 



point sources and a (sufficiently large) fixed physical aper- 
ture for extended sources. Since we expect a significant 
contribution from extended emission in the real data (e.g., 
Fig. 2), we focused on fixed physical apertures with proper 
radii of i? ap = 10, 20, 30, 40, and 50 kpc and background 
annuli extending over the larger of 30-60 kpc or 7? ap -2i? ap . 
Adopting a signal aperture of i? ap = 30 kpc (at the mean 
redshift of each bin) led to the most accurate reproduction 
of the input profiles with the highest S/N; our detection 
criterion is S/N > 1. This fixed physical aperture corre- 
sponds to angular apertures of ~ 15" — 5" over the range 
of mean redshifts we analyze (0.1 < z < 0.5, see Table 2 
and Fig. 2). The size of the field-averaged PSF begins to 
exceed that of significantly smaller apertures, particularly 
in the higher redshift bins (see Fig. 2). For point sources 
we could apply a redshift-dependent aperture correction 
to remedy this problem, but it is impossible to do so for 
extended sources when we are uncertain of their intrinsic 
profiles. For larger apertures, our signal-to-noise ratio be- 
gins to drop due to the increasing fraction of background 
photons included in the signal region. In the Monte Carlo 

30 



simulations, small deviations 10 kpc) from i? a 
kpc have little effect on the results or the signal-to-noise 
ratio but larger (factor of 2) changes in the aperture lead 
to serious problems. 

Using i? ap = 30 kpc and 30-60 kpc background an- 
nuli, which we adopt throughout the remainder of the pa- 
per, we successfully recover the input luminosity evolution 
(Eqn.[5]); i.e., when we fit our results with a double power 
law model (Eqn. [5]), we find scalings for the soft and hard 
band X-ray luminosities that are statistically consistent 
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Fig. 5. — The specific X-ray luminosity of all AGES sources 
compared to that of the normal galaxy sample. The solid line is a 
fit to the sum of the mean luminosity of the normal galaxies and 
the mean luminosity of the 58 X-ray detected AGN (see §2.2 and 
Tables 1 and 3) that lie within the absolute magnitude and redshift 
ranges we consider. The binned data shows the mean luminosity 
evolution of the normal galaxies alone; the point type conventions 
are the same as those used in Fig. \3\ The two points at each value 
of Lx demonstrate the evolution from the low to high redshift bin 
in each absolute magnitude strip (Fig. 1). 



with the input evolution rates and normalization: 



ix40,soft = (4.7 ±1.7) 



£ X 40,hard = (4.0 ± 2.4) 




1.0±0.2 



(1 + *) 



2.7±0.8 



0.9±0.3 



1 + 2 



,3.4±1.1 



(12) 



4. RESULTS 

We now turn to the real sources, starting with a com- 
parison of the X-ray emission from normal galaxies and 
from spectroscopic AGN. Fig. [5] shows the specific X-ray 
luminosity, (i X 4o) / {Lk /Lk*), as a function of (Lk/Lk*) 
for both the normal galaxies and the full sample includ- 
ing the X-ray and spectroscopically flagged AGN. With 
the top and right axes, we also show the inferred scaling 
with stellar mass, based on Eqn. ([3]), although this conver- 
sion does not apply when the bright AGN dominate the 
iiT-band flux. While the AGN have significantly larger X- 
ray luminosities than the "normal" galaxies, the "normal" 
galaxies are clearly detected. The two points at each Lk 
show the evolution of the normal galaxy X-ray emission 
from the low to high redshift bin in each absolute magni- 
tude strip (Fig. 1). When we fit the normal galaxy data 
as a double power law in K-band luminosity and redshift 
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Fig. 6. — Background-subtracted radial emission profiles of 360 
late-type (top) and 128 early-type (bottom) galaxies with z < 0.1, 
a mean redshift of (z) fs 0.07 and within 10' of CXO field centers. 
The solid, dashed, and dotted curves represent the total, hard, and 
soft band profiles, respectively, in each panel. The dot-dashed curve 
shows an empirical model of the PSF, based on the radial emission 
of « 500, z £ 1 AGN within 10' of a CXO field center. We have 
shifted the normalization of the PSF model to match the amplitude 
of the galaxy profile in each panel. 
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Fig. 7. — Redshift evolution of the Hardness Ratio (HR) for the 
late-type (squares) and early-type (circles) samples of galaxies. The 
late-type galaxies exhibit slightly harder spectra than the early-type 
galaxies, which is most likely due to the early-types' characteristi- 
cally larger halos of soft X-ray emitting gas. Typical HRs of HMXBs 
and LMXBs (dashed line) as well as Type I (dotted line) and Type 
II AGN (typically above the dashed line) have been included for 
reference. 



(Eqn. 03), we find that 

1*40 = (6.0 ± 1.0) 



£x40,hard = (4.6 ± 0.9) 



£x40,soft = (2.0 ± 0.5) 




l.lrbQ.2 



(1 + *) 



3.2±0.l 



1.2±0.3 



(1 + *) 



3.5±1.0 



1.1±0.2 



(1 + z) 2 ' 7±u - 7 (13) 



for the total, hard, and soft X-ray bands, respectively, in 
units of 10 40 ergs s -1 . Because our results indicate that 
L x is roughly proportional to Lk in all X-ray bands, the 
specific X-ray luminosities of the galaxies should evolve 
similarly regardless of their optical luminosity or mass. 

The X-ray luminosities we measure are roughly an or- 
der of magnitude brighter than the contribution expected 
from LMXBs (Eqn. [J), so we are left with three major 
candidates for the origin of the emission - HMXBs asso- 
ciated with recent star formation, accretion onto super- 
massive black holes (AGN), and hot gas. We use three 
tests to try to determine the dominant source(s). We 
first use fits to the AGES spectra to divide the normal 
galaxies into late-type and early-type subsamples. Doing 
so allows us to compare the X-ray emission from galax- 
ies with and without (optical) spectroscopic signatures of 
star formation. The resulting subsamples contain 3178 
late-type galaxies and 2968 early-type galaxies. Second, 
we examine the radial emission profiles of the low redshift 
galaxies, since emission due to star formation and hot gas 



will be more spatially extended than emission from AGN. 
Third, we compare the hardness ratios of the galaxies, in 
terms of the hard and soft band counts Chard and C so ft, 
HR = (Chard - C soft )/ (Chard + C^oft), to the typical values 
expected for X-ray binaries, AGN and hot gas. 

4.1. Radial Emission Profiles 

Fig. [HI shows the average, background-subtracted soft, 
hard, and total band emission profiles of 488 (360 late-type 
and 128 early-type) low redshift (z < 0.1), galaxies within 
10' of the CXO optical axis, excluding X-ray sources that 
violate our luminosity limit (Eqn. [HJ| and AGN identified 
through optical spectroscopy. The mean redshift of these 
galaxies is (z) ~ 0.07. Rather than use a model for the 
PSF (Eqn. [T0|) . we determined it empirically by averaging 
the emission profiles of w 500 z > 1 X-ray selected AGN 
from the XBootes survey. As with the galaxies, we used 
only AGN sources within 10' of the Chandra field centers 
and stacked the emission based on the optical positions. 
We then applied the angular corrections needed to put the 
PSF model on the same physical scale as the galaxies. 

We find that the emission profiles of the late-type galax- 
ies are clearly extended while those of the early-type galax- 
ies have a more pronounced nuclear component. We use 
the Kolmogorov-Smirnov (K-S) test to compare the ra- 
dial distribution of (total band) photons from the late- 
type galaxies, early-type galaxies, and high redshift AGN. 
Within 5 kpc of the stacking centers, we find that the 
early- type galaxy distribution (oc r~ 1A ) is consistent with 
the AGN distribution (K-S likelihood 97%), while the late- 
type galaxy distribution (oc r~ - 4 ) is not (K-S likelihood 
3%). However, if we analyze the full extraction region, 
neither galaxy profile is consistent with the AGN model 
(K-S likelihood < 1%). Between 10 and 30 kpc, both 
galaxy types show significant (^ 2.5 a) emission, which 
we attribute primarily to diffuse gas and LMXBs. Soft 
X-ray emitting gas has been observed out to > 20 kpc 
around late-type (e.g., Strickland et al. 2004) and early- 
type (e.g., Mathews & Brighenti 2003) galaxies, as have 
LMXBs (late-type: e.g., Wang et al. 1999; early-type: 
e.g., Kim et al. 2006), leading to extended profiles like 
those in Fig. 6. Our results imply that the emission from 
the low-redshift early-type galaxies comes from a combi- 
nation of AGN, hot gas, and LMXBs. In contrast, because 
the emission from the late- type galaxies is significantly less 
centrally concentrated than the emission of AGN - even 
within a radius of 5 kpc - the dominant source of the late- 
type galaxy flux is most likely a radially extended distri- 
bution of HMXBs and LMXBs. 

4.2. Hardness Ratios 

In Fig. [7] we show the hardness ratios (HRs) based on the 
hard and soft band counts of the two galaxy classes. To 
improve the statistical significance of our HRs, we combine 
hard and soft counts from bins of comparable rcdshifts in 
adjacent absolute magnitude strips (Fig. 1). For reference, 
we compare these HRs to those observed in these bands for 
accreting binaries (HR 0, e.g., Muno et al. 2004), Type 
I AGN (HR ~ -0.5, e.g., Rosati et al. 2003; Franceschini, 
et al. 2005), Type II AGN (0 < HR < 1, e.g., Rosati et 
al. 2003), and hot gas with a temperature of 0.5-1 keV 
(HR ~ -1). 
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In general, we observe that the late-type galaxies have 
slightly harder spectra than the early-type galaxies and 
that the spectra of both galaxy types become somewhat 
harder with increasing redshift. This is most likely a con- 
sequence of the larger amount of soft X-ray emitting gas in 
earlier type galaxies and the steadily increasing contribu- 
tions of HMXBs and AGN, respectively to the late-type 
and early-type galaxy spectra at higher redshifts. 

4.3. Luminosity Evolution 

To isolate the dependence of L x on Lk , we examined the 
evolution of the X-ray luminosity of the early and late- type 
galaxies in three narrow redshift slices: (z) = 0.20 ± 0.02, 
(z) = 0.25±0.02, and (z) = 0.35±0.03. As shown in Fig.[SJ 
the specific X-ray luminosities are relatively flat in Lk, 
indicating an approximately linear correlation between L x 
and Lk as found previously (e.g., David, Jones, & Forman 
1992; Shapley et al. 2001). If we fit these results with a 
power law model, L x oc L K , we find consistent values of 
a among the three redshift slices for both the late-type 
( a htrd E = 1-19±0.31; a^4 TE = 1-12±0.15) and early-type 
(«wd LY = 1-09 ± 0.23; af f t RLY = 1.07 ± 0.15) galaxies. 

Figs. 191 and ITTJ1 illustrate the dependence of the specific 
X-ray luminosity of the late-type and early-type galaxies, 
respectively, on redshift (left panels) and optical luminos- 
ity (right panels). If we fit the late- type galaxies with the 
double power-law model (Eqn. [5]), we find that, in units of 
10 40 ergs s -1 , 

iZfflLi = (5-9 ± 1.6) (^-) 1 ' 2± °' 4 (1 + z) 3 - 6±1 - 2 , 

L x L 4 A ™ ft = (2.5 ± 0.6) 1 ' 1± °' 2 (1 + z)«±"<14) 

The reduced x 2 (per degree of freedom) for these fits are 
0.78 and 0.81 for the hard and soft bands, respectively. If 
we fit the early- type galaxies, we find that 

/j x 1.0±0.3 

= (3-0 ± 0.6) (gl J (1 + z) 3 ™, 

Llt^l = (1.8 ± 0.4) (|*L) 11±0 ' 2 (1 + z) 31±0 8 (15) 

In this case, the reduced % 2 values are 0.86 and 1.14 for the 
hard and soft bands, respectively. The reduced % 2 values 
for the best-fit curves given in the remainder of the paper 
are all similarly good, i.e., « 1 ± 0.2. In all cases, the 
upper limits set by non-detections (bins with S/N < 1) 
are used in conjunction with detections to establish the 
best-fit curves. In particular, if any best-fit curves exceed 
these upper limits, the squared difference is added to the 
X 2 sum. 

The results for the individual luminosity and redshift 
bins arc summarized in Tables 2 and 3. In spite of the evi- 
dence that different physical processes dominate the emis- 
sion of the late-type and early-type galaxies, the two pop- 
ulations exhibit qualitatively similar behaviors. In partic- 
ular, 

• the X-ray luminosity monotonically increases from low 
mass galaxies at low redshift to high mass galaxies at high 
redshift and 

• the X-ray luminosity is significantly higher than what 
would be expected from LMXBs on the basis of the K- 



band luminosities, Eqn. @, assuming a hardness ratio of 
HR = (e.g., Muno et al. 2004). 

These similarities are not surprising, as we expect the 
X-ray emission of LMXBs and hot gas to be nearly time- 
independent for z <, 0.5, while the HMXB and AGN emis- 
sion should be evolving rapidly and at approximately the 
same rate (oc (1 + z) 3±1 , e.g., David, Jones, & Forman 
1992; Shapley et al. 2001, Hogg 2001, Norman et al. 2004, 
and Schiminovich et al. 2005 and BaOl and Ba05, respec- 
tively) . 

4.4. Comparisons to Previous Work 

In general, we find excellent agreement between our 
results and those of earlier studies of normal galaxy X- 
ray evolution, but our uncertainties are smaller. There 
are three general categories of earlier results: surveys of 
nearby individual galaxies, stacking analyses of low to 
intermediate-redshift galaxies and stacking analyses of in- 
termediate to high-redshift galaxies in small-area, deep 
fields (CDF-N and CDF-S). Taking the late-type and early- 
type galaxies in turn, we first consider the Shapley et al. 
(2001, S01) and O'Sullivan et al. (2001, OS01) surveys of 
local spiral and early- type galaxies. Second, we compare 
our results to the Georgakakis et al. (2003, G03) stacking 
analysis of 2dF galaxies at z ~ 0.1 and to our own earlier 
stacking analysis of massive early-type galaxies at 0.3 <, 
z <y 0.9 (B05). Finally, we consider the Hornschemeier et 
al. (2002, H02) and Lehmer et al. (2007) stacking analy- 
ses of normal galaxies in the CDF-N at 0.4 <J z <^ 1.5 and 
optically bright early-type galaxies in the extended CDF-S 
out to z 0.7, respectively. 

Shapley et al. (2001, S01) studied the correlation be- 
tween the X-ray (0.2-4.0 keV) and ff-band flux of 234 
local spiral galaxies with no sign of AGN activity. Their 
results indicate a trend of 

^s E oi(- * 0) * (2 ± 1) 1 ' 2± °' 1 , (16) 

assuming H— K= 0.27 mag and a 50% uncertainty in the 
normalization of the correlation. It is clear that Eqn. (jTBJ) 
agrees very well with Eqn (|14|) at z = 0, and, as shown in 
Fig. [51 the agreement extends to all redshifts if we allow 
for evolution (see Fig. 9). Similar trends have been found 
in other local surveys (e.g., David, Jones, & Forman 1992). 

Georgakakis et al. (2003, G03) considered 200 galaxies in 
the 2dF redshift survey with (Lg) ~ at a mean redshift 
(z) ~ 0.1. They measured mean soft band luminosities of 
log(L x /ergs s" 1 ) ~ 40.2 ± 0.2, 39.7 ± 0.4 and 40.5 ± 0.2 for 
the combined, late-type and early-type samples. They also 
estimated an upper limit of log(L x /ergs s _1 ) SS 40.5 on the 
hard band luminosities of the samples. These results are 
consistent with ours, particularly if we allow for positive 
redshift evolution of the X-ray luminosity between (z) ~ 
0.1 for G03 and our points for galaxies of similar optical 
luminosities at z ~ 0.25 and 0.45. 

Hornschemeier et al. (2002, H02) conducted a stack- 
ing analysis of late-type galaxies in the HDF-N/CDF-N, 
finding that 29, ~ 0.7L, galaxies in the redshift range 
0.4 — 0.75 have mean total (0.5 — 8 keV) and soft band 
(0.5-2.0 keV) luminosities of L x40 = 2.9±0.4and 1.3±0.2, 
respectively. These too are consistent with the values we 
find for ~ 0.7L* galaxies at 0.2 ^ z ^ 0.3 (3.6 ± 0.8 and 
1.6 ± 0.3 respectively). 
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Fig. 8. — The specific X-ray luminosity evolution of the early and late- type galaxies in three narrow redshift slices: (2) = 0.20 ± 0.02 
(three-pointed stars), (2) = 0.25 ± 0.02 (five-pointed stars), and (2) = 0.35 ± 0.03 (seven-pointed stars). In all redshift slices, the specific 
X-ray luminosities are relatively flat in Lk (roughly, L x ,latc oc L^ 2±0 ' 3 oc L x ,carlyi see §4.3), indicating an approximately linear correlation 
between L x and Lk as found previously (e.g., David, Jones, &: Forman 1992, Shaplcy ct al. 2001). 
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Fig. 9. — Left: The mean hard and soft band specific X-ray luminosities of the late-type sample of galaxies as a function of redshift. The 
point type conventions are the same as those used in Figure|3] The best-fit curves for the hard and soft band emission are also shown. Right: 
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scatter in and error bars on the S01 data points. 
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Fig. 10. — Same as Figure [9] for the early-type sample of galaxies. For comparison, we have plotted the specific X-ray luminosities of 
2 ~ 0.1 L* galaxies found by Georgakakis et al. (2003, G03, 4-pointed star in the soft band panel; upper limit in the hard band panel) and 
Brand et al. (2005, B05, 6-pointed star in both panels at z = 0.4 ± 0.1). The best-fit curve for soft and hard band emission from the AGES 
galaxies is also shown. The shaded region in the lower right panel shows the local O'Sullivan et al. (2001, OS01) L x soft^-K") trend, roughly 
log(L x 4o) = (0.3 ± 0.3) + (1.7 ± 0.2) log(L k / L k t) , which encompasses both the scatter in and error bars on the OS01 data points. 
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Table 1 

X-ray Counts From Stacked Normal Galaxies & Excluded AGN 



Class 


Source 


Number of 


Net Hard 


Net Soft 






Objects 


Counts 


counts 


Normal Galaxies 


Late- Type Galaxies 


3178 


240 


327 




Early- Type Galaxies 


2968 


246 


342 


Spectroscopic AGN 


Broad Line AGN 


27 


319 


872 




Narrow Line AGN 


78 


496 


404 


X-ray Sources 


Broad Line AGN 


2d 


319 


872 




Narrow Line AGN 


34 


468 


367 
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Table 2 

X-ray Luminosity Upper Limits for the Optically Faintest AGES Galaxies 







(z) -^Vg a i ax i cs 


L X 40 


(mag) 


(mag) 


(h/s) 



Late- Type Galaxies 



-17.5 


-18.5 


0.0693 


126 


< 0.6/ 


< 0.2 


-17.5 


-18.5 


0.0951 


126 


< 1.9/ 


< 1.0 


-18.0 


-19.0 


0.0969 


228 


< 1.2/ 


< 0.7 


-18.0 


-19.0 


0.1286 


228 


< 1.2/ 


< 0.4 




Early- Type Galaxies 




-17.5 


-18.5 


0.0756 


13 


< 1.8/ 


< 0.4 


-17.5 


-18.5 


0.0882 


13 


< 2.2/ 


< 0.9 


-18.0 


-19.0 


0.1051 


59 


< 2.6/ 


< 0.6 


-18.0 


-19.0 


0.1303 


56 


< 3.4/ 


< 0.9 


-18.5 


-19.5 


0.1166 


119 


< 2.4/ 


< 0.4 


-18.5 


-19.5 


0.1586 


119 


< 3.5/ 


< 1.5 



Note. — The range of absolute ij-band magnitudes of the 
galaxies in a bin is M^ max > Mr > Afp,. m i n and the mean 
rcdshift is (z). -/V ga i ax ; os is the effective number of galaxies 
in the bin, including the correction for sparse sampling. The 

1 a upper limits on the hard (h; 2—7 keV) and soft (s; 0.5— 

2 keV) band X-ray luminosities, L x 4o, of the optically faintest 
AGES galaxies are given in units of 10 40 ergs s _1 in the final 
column. 



Table 3 
Summary of Results 



(M R )/(M K ) 


<*> 


-^galaxies 


Net Counts 


S/N 


ix40 


(L x 4o)/{Lk /Lk*) 


(mag) 




(h/s/t) 


(h/s/t) 


(h/s/t) 


(h/s/t) 








Late-Type 


Galaxies 






-18.98/-21.31 


0.1159 


375 


17/ 34/ 51 


0.8/2.0/1.8 


0.7/ 0.4/ 1.1 


4.3/ 2.5/ 6.8 


-19.02/-21.33 


0.1614 


376 


- 1/ 17/ 16 


— /1.3/0.7 


< 1.4/ 0.4/ 0.3 


< 9.2/ 2.5/ 1.8 


-19.42/-21.75 


0.1434 


580 


-4/ 71/ 67 


— /3. 9/2.3 


< 0.8/ 0.7/ 0.6 


< 3.9/ 3.4/ 2.8 


-19.47/-21.80 


0.1979 


576 


27/ 52/ 79 


1.5/3.6/3.4 


2.1/ 1.1/ 3.2 


9.4/ 5.1/14.5 


-19.93/-22.25 


0.1709 


667 


55/ 61/116 


2.4/3.3/4.0 


2.6/ 0.8/ 3.5 


7.7/ 2.4/10.1 


-19.94/-22.26 


0.2464 


667 


- 23/ 45/ 22 


— /3. 4/1.0 


< 2.6/ 1.4/< 3.0 


< 7.5/ 4.0/< 9.0 


-20.36/-22.68 


0.2032 


555 


39/ 59/ 98 


2.1/4.2/4.2 


3.4/ 1.4/ 4.8 


6.3/ 2.6/ 8.9 


-20.44/-22.75 


0.3011 


556 


3/ 30/ 33 


0.2/2.7/1.9 


0.6/ 1.8/ 2.4 


1.1/ 3.3/ 4.4 


-20.85/-23.16 


0.2465 


403 


46/ 50/ 96 


3.1/4.3/5.1 


8.5/ 2.5/11.0 


10.2/ 3.1/13.2 


-20.90/-23.20 


0.3774 


403 


36/ 34/ 70 


3.6/3.7/5.1 


17.5/ 4.6/22.2 


21.1/ 5.6/26.6 


-21.27/-23.57 


0.2878 


195 


25/ 39/ 64 


2.4/4.4/4.7 


13.2/ 5.9/19.1 


10.2/ 4.5/14.8 








Early- Type 


Galaxies 






-19.55/-21.88 


0.1482 


328 


14/ 31/ 45 


0.8/2.2/2.0 


1.0/ 0.6/ 1.6 


4.5/ 2.8/ 7.3 


-19.51/-21.84 


0.2047 


327 


11/ 19/ 30 


0.7/1.7/1.6 


1.6/ 0.8/ 2.4 


7.3/ 3.6/10.9 


-19.99/-22.31 


0.1672 


560 


47/ 62/109 


2.1/3.4/3.8 


2.6/ 1.0/ 3.6 


7.6/ 2.8/10.4 


-20.00/-22.32 


0.2438 


561 


20/ 53/ 73 


1.2/3.9/3.4 


2.5/ 1.9/ 4.4 


7.4/ 5.5/12.8 


-20.45/-22.77 


0.2102 


804 


31/ 72/103 


1.3/4.0/3.5 


2.0/ 1.3/ 3.3 


3.7/ 2.4/ 6.1 


-20.45/-22.76 


0.3146 


796 


37/ 65/102 


2.1/4.5/4.5 


6.0/ 3.0/ 8.9 


11.2/ 5.5/16.7 


-20.90/-23.21 


0.2564 


844 


38/103/141 


2.0/6.1/5.6 


3.6/ 2.8/ 6.4 


4.4/ 3.3/ 7.7 


-20.90/-23.20 


0.3823 


844 


55/ 98/153 


3.4/6.7/7.0 


13.3/ 6.6/19.9 


15.9/ 8.0/23.9 


-21.33/-23.63 


0.2954 


594 


45/ 69/114 


3.0/4.9/5.5 


8.4/ 3.6/12.0 


6.5/ 2.8/ 9.2 


-21.38/-23.67 


0.4691 


593 


48/ 42/ 90 


3.8/3.7/5.3 


26.5/ 6.6/33.1 


20.4/ 5.1/25.5 



Note. - - The mean i?-band and ii"-band absolute magnitudes of the galaxies in a bin are (Mr) and (Mk) 
and the mean redshift is (z). iVg a i ax [ es is the effective number of galaxies in the bin, including the correction for 
sparse sampling. The Net Counts column gives the number of photons observed above background in the hard (h; 
2—7 keV), soft (s; 0.5—2 keV) and total (t; 0.5—7 keV) bands; the corresponding signal-to-noise ratio (S/N) of the 
measurements appears in the adjacent column. (Our detection criterion is S/N > 1). L x 4o is the mean X-ray 
luminosity in the hard, soft and total bands in units of 10 40 ergs s — 1 , and the final column is the specific X-ray 
luminosity, (£ x 4o,(h,s,t) )/{I j K /Lk *), for each band. In the cases of negative net counts, 1 a upper limits have been 
provided. Only bins that contain > 100 galaxies and have net stacked emission with S/N > 1 in at least one X-ray 
band have been included in the table. 
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Fig. 11. — The mean hard (top) and soft (bottom) specific X-ray 
luminosities of ~ L* early-type galaxies from G03 (2dF, z = 0.1), 
this work (AGES: filled points), and B05 (NDWFS: open points at 
z <; 0.4). The solid lines are global fits to the combined samples 
(see §4.4). 



Turning to the early-type galaxies, we first checked the 
relative contributions of hot gas and AGN emission by 
comparing our results to the ROSAT PSPC (0.1-2.5 keV) 
X-ray survey of 401 local early- type galaxies discussed 
in O'Sullivan et al. (2001). Because the hot gas emis- 
sion from the early-type galaxies is primarily a function of 
galaxy mass rather than epoch, there should be little dif- 
ference between the soft band flux of the local and AGES 
galaxies (at 0.15 <; (z) <; 0.45), provided that the soft 
band AGN contamination in the AGES sample is small. 
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Fig. 12.— The Specific Star Formation Rate (SSFR = 
(SFR)/(M*)) as a function of (Lk) and (M«) based on the hard 
(dashed) and soft (dotted) X-ray luminosities of the AGES late- 
type galaxies and Eqns. JTJ - For the soft band values, the 
point type conventions are the same as those used in Fig. 3. We 
compare our results to those of Bauer et al. (2005) (six-pointed 
stars), and Bell et al. (2005) (four-pointed star). 
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Fig. 13. — Estimates of the star formation rate density p* based on 
the hard (solid) and soft (dashed) X-ray luminosities of the AGES 
late-type galaxies as compared to the GALEX results (Schiminovich 
et al. 2004) after either their minimum (A^ v ) or best-fit (A™{^) 
corrections for dust extinction have been applied. 
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Fig. 14. — The specific supermassive black hole (SMBH) accretion 
rate as a function of redshift based on the hard band X-ray lumi- 
nosity of AGES early-type galaxies (in excess of LMXB emission) 
and Eqn. Q. The dashed line is the curve of best fit. The point 
type conventions are the same as those used in Fig. 3. 

Fig.fTOlshows that the soft band luminosity from the AGES 
early-type galaxies is, in fact, comparable to that of the lo- 
cal sample even at high optical luminosities (which tend to 
be associated with higher redshifts in the AGES sample). 
In contrast, the total luminosity of the AGES early- type 
galaxies is consistently higher than that of the local sam- 
ple. Fig.[Tn]also demonstrates that the hard X-ray compo- 
nent of the early-type galaxies is well in excess of the esti- 
mated LMXB emission and undergoes rapid, positive red- 
shift evolution. These findings reinforce our conclusions in 
§4.2, that the soft band flux of the AGES early-type galax- 
ies comes primarily from hot gas, while the more rapidly 
evolving hard component comes primarily from AGN. 

Our results for the massive (~ L*) early- type galaxies 
agree with those of G03 and with our earlier work (B05), 
where we found that bright red (Mr < —21.3) galaxies at 
z ~ 0.4 had total X-ray luminosities of L x ~ 10 41 ergs s . 
By incorporating the results of G03 and B05 with ours, we 
gain more leverage on the behavior of L* early-type galax- 
ies at lower (z ~ 0.1) and higher (0.5 <S z <^ 0.9) redshifts 
than what we can probe with AGES alone. Using this 
combined sample, we find that the hard and soft band 
luminosities of L» early-type galaxies follow the trends 

^ffid(£.) = (3-6 ± 1.5)(1 + zf- 5±1 - 7 (17) 

and 

£x 4 A o R soft(£*) = (2.0 ± 0.8)(1 + z) 2 -^ 1 ' 1 (18) 

out to z ~ 1. These global fits, which are shown in Fig.QTJ 
are consistent with those found for the AGES L* early- 
type galaxies alone (Eqn. [TBI at Lk — Lk*)- 




^0 0.1 0.2 0.3 0.4 0.5 

Z 

Fig. 15. — The comoving supermassive black hole (SMBH) accre- 
tion rate density as a function of redshift for the AGES early-type 
galaxies. Because we are dealing with normal galaxies rather than 
X-ray bright AGN, our results lie roughly a factor of four below the 
Ba05 data (open squares at 0.15 < z < 0.3 and 0.35 < z < 0.5) and 
the (solid) Ba05 best-fit curve (which is based on additional data 
out to z ~ 1.2). The slopes of the best-fit curves shown do agree to 
within errors, however. Our value is 2.9 rfc 0.7 (dashed line), while 
Ba05 finds 3.2 ± 0.8 (for z < 1.2). 



Lehmer et al. (2007) conducted a stacking analysis of 
222 optically bright (~ L#*) early-type galaxies in the E- 
CDF-S with photometric redshifts ranging from to 0.7. 
The mean hard X-ray luminosity trend inferred for AGN in 
the Lehmer et al. galaxy sample is in excellent agreement 
with the hard band trend found in B05 (see Lehmer et al. 
2007; Fig. 13) and in the present paper (Fig.lPPj). 

5. INTERPRETATION AS STAR FORMATION AND NUCLEAR 
ACCRETION 

In this section we interpret the trends in X-ray luminos- 
ity in terms of star formation (HMXBs) and nuclear ac- 
cretion (AGN) after estimating and subtracting the addi- 
tional contributions of LMXBs and hot gas. Under the as- 
sumption, discussed in the introduction, that Eqn. [2] holds 
out to z ~ 0.5, LMXBs represent a modest contribution 
to the emission from the late-type (sS 15%) and early-type 
(<J 20%) galaxies. The comparison of the (primarily) hot 
gas emission from local early-type galaxies (O'Sullivan et 
al. 2001) shown in Fig. [TOl suggests that, at all redshifts, 
~ 50% of the soft band (and <, 20% of the total) emis- 
sion of the AGES early-type galaxies is produced by hot 
gas. Based on these arguments, we conclude that <; 85% 
of the X-ray flux of the late-type galaxies is produced by 
HMXBs and £ 60% of the total (and £ 80% of the hard) 
X-ray flux of the early-type galaxies is produced by low 
luminosity AGN. In the analysis that follows, we use the 
hard and soft band luminosities of the late-type galaxies 
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to estimate the evolution of star formation rates in nor- 
mal galaxies and the hard band luminosity of the early- 
type galaxies to estimate the growth rates of supermassive 
black holes in normal galaxies. 

Before we can implement these assumptions, we must 
also consider the possibility that the flux from the high 
redshift late-type galaxies is being contaminated by AGN. 
Since we cannot resolve their radial emission profiles as we 
did for the nearby subsamples (§4.1) and their luminosity 
evolution is fairly similar to that of the early-type galax- 
ies, which we largely attribute to AGN, we must find other 
ways to address this possibility. Fortunately, there are 
three key properties of the late-type galaxy emission that 
suggest the level of AGN contamination is small. First, 
the star formation signatures in the late-type galaxy op- 
tical spectra are necessarily correlated with a significant 
flux from HMXBs (Eqn. 1), suggesting that only a small 
fraction of the late-type galaxy emission could come from 
AGN. Second, the X-ray emission from late-type galaxies 
is extended at low redshifts (Fig. 6), directly demonstrat- 
ing the absence of AGN among the nearby members of 
the sample. Third, the Poisson and bootstrap error esti- 
mates are similar at all redshifts, which they would not be 
if a small number of individual, high-count sources (AGN) 
were hidden within the sample. This disfavors small to 
moderate AGN contamination. Significant AGN contami- 
nation is also disfavored by our analysis and comparisons 
to previous work in §5.1. 

5.1. Star Formation Rates in the Late- Type Galaxies 

Studies of star formation over the redshift range <S z 
3 have demonstrated that the relationship bewteen X-ray 
luminosity and star formation rates is reasonably univer- 
sal (e.g., Brandt et al. 2001; Seibert et al. 2002; Nandra 
et al. 2002; Cohen 2003; Grimm et al. 2003; Persic et al. 
2004; Reddy & Steidel 2004; Lehmer et al. 2005). We use 
Eqn. dJ) to estimate the hard band SFR from = 
^Ifd - ffl, where L™ is based on Eqn . ( 2) and 
a typical LMXB spectrum with H R = (e.g., Muno et 



al. 2004) and L^Ird is S iven in Ecm - 



SFR 



lM , d = (6A±l.8)(L K /L K *) 1 - 2±0A 
x(l + z) 3 - 4±1 - 2 M yr- 1 . 



3). The result is 



(19) 



Using Eqn. © to estimate the stellar mass M* of the 
galaxies, we find that SFRhard corresponds to a specific 
star formation rate (SSFR = SFR/ A/*) of 



SSFRhard — 



14 ±2.4) x 10~ 2 (1 + z) 3 ' 4±1 ' 2 
x(M,/10 11 M o ) n - 2±a4 Gyr- 1 . 



(20) 



Based on the work of Ranalli et al. (2003, as modi- 
fied by Gilfanov et al. 2004) we also use Eqn. |T]) to es- 
timate the soft band SFR from L™xb = L late_ l lmxb : 



SFR soft = (3A±0.7)(L K /L K , 



il.l±0.2 



x(l 



\2.4±1.0 



Meyr" 1 



(21) 



and a corresponding specific star formation rate of 

SSFR soft = (4.3 ± 0.9) x 10~ 2 (1 + z) 2A±1M 

x(Af,/lO 11 M )°- 1±a2 Gyr- 1 . (22) 
As shown in Fig. [TJl the results of earlier SSFR studies 
by Bauer et al. (2005), Bell et al. (2005), and Feulner et 
al. (2004) fall within the range encompassed by our hard 
and soft band SSFR trends fEqns.H^TJl and |2"2"| . but we have 
smaller uncertainties due to our larger sample size. Our 
results also agree with those of more recent analyses of 
comparably large galaxy samples (Noeske et al. 2007 and 
Zheng et al. 2007) in bins of similar stellar mass (10.5 <, 
log(M*/M Q ) & 11) and redshift (0.2 <, z <, 0.45 and z ~ 
0.3, respectively). 

We compute the star formation rate per comoving vol- 
ume, p*, by integrating Eqns. (fT§|) and (f2"Tj) over the co- 
moving density of the late-type galaxies (see §2.2). These 
results are shown in Fig. [T3] We find 



P*,hard = (2.1 ±0.4) x 10- 2 (l + z) 
and 

p*, Boft = (1.6 ±0.3) x 10" 2 (1 



2.9±0.7 



Mr. 



yr Mpc 3 

2.5±0.5 M © 



(23) 



(24) 



yr Mpc 3 

for the hard and soft bands respectively. As illustrated in 
Fig. 1131 our results are broadly consistent with the range of 

values inferred from the standard extinction-corrected 
ultraviolet luminosity measurements of GALEX (Schimi- 
novich et al. 2004). 

It is important to note that significantly larger extinc- 
tion corrections have already been ruled out by neutrino 
detectors. Values of p* more than 50% larger than p*,hard 
would be in serious conflict with the Super-Kamiokande 
upper limit (Malek et al. 2003) on the Diffuse Supernova 
Neutrino Background (DSNB) from M > 8M© stars (Stri- 
gari et al. 2005, Hopkins and Beacom 2006). 

Our soft band estimates of p* sit closer to the GALEX 
minimum extinction correction curve. Nevertheless, the 
consistency of our hard and soft band results with each 
other and with the range of GALEX results indicates that 
we are in general agreement with earlier studies (e.g., based 
on Ha measurements: Tresse & Maddox 1998, Tresse et 
al. 2002, Perez- Gonzalez et al. 2003, and Brinchmann et 
al. 2004 and the meta-analysis of Hogg 2001) as well as 
more recent comprehensive studies (e.g., Hopkins & Bea- 
com 2006). 

The level of agreement between our hard and soft band 
results also provides a check on the contamination of our 
sample by AGN. We first note that the soft band SSFR and 
p» values (Eqns. [HI and [M]) are systematically lower than 
the corresponding hard band values fEqns.l2"0"landl2"3"]). but 
this is not uncommon. Many previous studies have also 
generally inferred lower SFR and p* values based on soft 
band flux (e.g., G03 and Norman et al. 2004), suggesting 
that a low-normalization conversion between the soft band 
X-ray flux and the star formation rate (e.g., Hornschemeier 
et al. 2005), rather than contamination in the hard band, 
is to blame for this discrepancy. Furthermore, if sources 
harder than X-ray binaries, like Type II AGN (see Fig. [7]), 
were contributing excess hard band emission to the late- 
type galaxy flux, p* ; hard (Eqn. |2"3"]) would lie above rather 
than slightly below the GALEX best-fit curve (Fig. [T3|). 
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5.2. Accretion Rates in the Early-Type Galaxies 

For the early-type galaxies, we use Eqn. © to estimate 
and subtract the contamination from LMXBs and use the 
remaining hard X-ray luminosity to determine the average 
supermassive black hole growth rate via Eqn. ^ . Assum- 
ing a constant accretion efficiency of e = 0.1, we find that 

A/bh - (5.8 ± 1.1) x lQ- & {L K /L Kjf ) l - x±0A (25) 
x(l + z) 3 - 3±a7 M yr- 1 . 

This growth rate (see Fig. ll4p is roughly an order of magni- 
tude below the estimates of BaOl (Mgg 1 £ lO^Mgyr" 1 ), 
because we are considering very different source popula- 
tions. The BaOl estimates are based on directly detected 
X-ray sources with luminosities exceeding 10 42 ergs s _1 , 
while we have eliminated such bright sources and consider 
only "normal" galaxies with mean luminosities of order 
10 41 ergs s . The rapid redshift evolution of the accre- 
tion rate in normal galaxies is consistent, however, with 
the redshift evolution of bright AGN found by BaOl and 
Ba05. 

We estimate the accretion rate density /?bh by inte- 
grating the accretion rate per unit (stellar) luminosity 
over the early-type galaxy luminosity function (see §2.2). 
We find that 

Pbh = (2.0±0.4)xl0" 8 (l+z) 2 ' 9±a7 M Q yr- 1 Mpc" 3 . (26) 

As shown in Fig. 1151 the evolution rate we determine is 
consistent with the rates found by BaOl (pbh oc (l + z) 3 ) 
and Ba05 (pbh oc (1 + z ) 3 - 2±0 - 8 ) but we have a systemati- 
cally lower normalization for pbh because of the different 
source populations we have examined. This suggests that 
even at z <J 0.5 the aggregate growth rate of black holes in 
normal galaxies is still a factor of a few less than that oc- 
curring in the much smaller number of bright, low-redshift 
AGN. 

6. SUMMARY AND CONCLUSIONS 

After eliminating AGN flagged by AGES spectroscopy 
and X-ray sources brighter than our luminosity limit (Eqn. 
8), we used a stacking technique to determine the av- 
erage X-ray properties of a magnitude-limited sample of 
6146 (R < 20 mag, z < 0.6) apparently normal galaxies. 
The results substantially improve our knowledge of normal 
galaxy X-ray evolution over the redshift range 0.1 ^ z <J 
0.5, which has been difficult to probe based on either local 
studies of individual objects or small-area deep fields. 

We spectroscopically divided the galaxies into late-type 
and early-type subsamples. The spectroscopic signatures 
of star formation, spatially extended radial emission pro- 
files, luminosity evolution, and low LMXB emission (<, 
15%) of the late-type galaxies suggest that their X-ray 
emission is dominated by HMXBs and therefore traces the 
star formation rate. 

Conversely, because the early-type galaxies lack spectro- 
scopic evidence for star formation, yet have rapidly evolv- 
ing hard X-ray luminosities that are well in excess of that 
expected from LMXBs (<, 20%), we concluded that their 
emission is increasingly dominated by AGN at higher red- 
shifts. 



When we use a double power law (L x oc L^(l + z) /3 ', 
Eqn. to fit the trends in the X-ray emission, we see 
that the X-ray luminosity increases monotonically from 
low mass galaxies at low redshift to high mass galaxies 
at high redshift for both galaxy types. The redshift evo- 
lution of the late-type galaxies, L x oc (1 + z) 3±1 , is in 
good agreement with previous estimates (e.g., H02, Ptak 
et al. 2001) and with theoretical expectations for normal, 
star-forming galaxies (Ghosh & White 2001). The optical 
luminosity dependence of the late-type galaxy emission, 
roughly L x oc L^ 2±0 ' 3 , also matches previous results (e.g., 
S01; David, Jones, & Forman 1992). In addition, the spe- 
cific star formation rates, SSFR = SFR/M*, and star for- 
mation rate densities, p*, we infer from the hard and soft 
band emission of the late-type galaxies span the range of 
previously reported values fFigs . [T2l and [T5|) . 

The specific X-ray luminosities we found for the early- 
type galaxies are in good agreement with the results of 
past studies (O'Sullivan et al. 2001, G03, B05, and Lehmer 
2007). Additionally, our work provides the first X-ray es- 
timates of the SMBH accretion rates, Mbh, and accretion 
rate densities, pbh, in normal galaxies at these redshifts. 
Our findings suggest that the redshift evolution of low lu- 
minosity AGN in our early-type galaxy sample is similar to 
that of higher luminosity AGN (BaOl and Ba05), which we 
exclude from our analysis. However, the lower luminosity 
AGN contribute little to the overall growth of supermas- 
sive black holes. In general, our work shows that there is 
a continuum rather than a sudden break in the star for- 
mation and SMBH accretion histories of galaxies from the 
powerful starbursts and AGN of the past to the fainter, 
optically-normal galaxies more prevalent today. 

Our analysis of the X-ray evolution of galaxies in the 
NDWFS XBootcs field can be significantly expanded in the 
future. A second phase of AGES has doubled the size of 
the spectroscopic sample and extended the redshift range 
of galaxies. If combined with photometric redshifts based 
on the very extensive, multiwavelength NDWFS photom- 
etry, it will be possible to reach z ~ 1 to a uniform lumi- 
nosity limit and greatly reduce statistical errors. 
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